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Abstract 

A previously developed mathematical model is amended to more accurately incorporate the effects of lift 
and drag on single dust particles in order to predict their behavior in the wake of high velocity gas flow. 

The model utilizes output from a CFD or DSMC simulation of exhaust from a rocket nozzle hot gas jet. An 
extension of the Safftnan equation for lift based on the research of McLaughlin (1991) and Mei (1992) is 
used, while an equation for the Magnus force modeled after the work of Oesterle (1994) and Tsuji et al 
(1985) is applied. A relationship for drag utilizing a particle shape factor ((f) = 0.8) is taken from the work 
of Haider and Levenspiel (1989) for application to non-spherical particle dynamics. The drag equation is 
further adjusted to account for rarefaction and compressibility effects in rarefied and high Mach number 
flows according to the work of Davies (1945) and Loth (2007) respectively. Simulations using a more 
accurate model with the correction factor (€ = 0.8 in a 20% particle concentration gas flow) given by 
Richardson and Zaki (1954) and Rowe (1961) show that particles have lower ejection angles than those that 
were previously calculated. This is more prevalent in smaller particles, which are shown through velocity 
and trajectory comparison to be more influenced by the flow of the surrounding gas. It is shown that 
particles are more affected by minor changes to drag forces than larger adjustments to lift forces, 
demanding a closer analysis of the shape and behavior of lunar dust particles and the composition of the 
surrounding gas flow. 
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Key Point #1 



Key Point #2 

Assumption 

Ballistics calculations ignore the 
temporal variation of the gas. 
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golf hall rolling across a p 
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Trajectory Models Simulated 



2 nd Order Taylor’s Series 
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Simple Model 

Vertical Gradient of Horizontal Force * 0 
=> Bernoulli Lift Force 
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CFD Results (Case C2) 

CFD Gas Velocity [m/s] CFD Gas Density [kg/m*] 
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Trajectory Dependence on 
Surface Topography 
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Estimated Dust Ejection Speed and 
Angle from Ballistics Simulations 
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Particle trajectory angle* 
relative to ground for various 
particle tire* and CFD cate* 


CFD model boundary 


Summary 

Point *1 


Angles Trajectory Angles 

of simulated derived from 

particle Apollo video 

trajectories analysis 


• Craters can generate a substantial disruption to 
particle ballistics motion 

• Particles can easily (all into craters. 

• Particle motion down-stream from craters, at the 
surface layer, may be suppressed 


Points #3 

Comparison of Integration Algorithms 

• I" 4 Order Taylor Series algorithm is very su 
and well suited for single particle simulations, m 
long as the step-aize is sufficiently small 

• Kerman's Method is Much Better 

• KM is superior to both 
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